1D-FFLO state in the absence of time-reversal symmetry breaking 
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A novel route to a one-dimensional Fulde-Ferrell-Larkin-Ovchinnikov (1D-FFLO) state in the 
absence of broken time-reversal symmetry is proposed in this paper. At present such a state may be 
encouraged in a clean AlAs quantum wire. Using the AlAs quantum wire as an example it is shown 
using bosonization and the renormalization group approach that the 1D-FFL0 state can arise due 
to a combination of Coulomb interactions and the unique bandstructure arrangement of the AlAs 
quantum wire. The present theoretical proposal is very general and is applicable to other systems 
with similar fermionic interaction terms. 

PACS numbers: 73.21.Hb,74.10.+v,74.20.-z, 74.20. Mn,74.78.Na 



A superconducting state with a finite pairing momen- 
tum had been long proposed by Fulde and Ferrell 
and Larkin and Ovchinnikov 0]. Presently this state is 
known as the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
(or LOFF) state. The FFLO state is widely realized from 
condensed matter systems to high energy physics [3|. In 
a FFLO state Cooper pairs with a nonzero center-of-mass 
momentum are encouraged in the presence of an exter- 
nal magnetic field. The central physical ingredient of the 
originally studied FFLO state is based on the idea that 
Cooper pairs in a spin-singlet superconductor are com- 
posed of fermions with opposite spins. In the presence 
of a magnetic field these electrons couple to the external 
applied field through the Zeeman coupling. This in turn 
tends to polarize the electrons along the direction of the 
magnetic field and allows the system to gain polarization 
energy. However, the pairing of opposite spins is favor- 
able for condensation energy. Due to these competing 
tendencies the superconducting state undergoes a transi- 
tion to the FFLO state with a finite pairing momentum 
and eventually enters the normal state as the magnetic 
field is further increased [l], Q • 

There have been several theoretical proposals that such 
a state may be achievable in heavy fermion systems and 
org anic and cuprate superconductors 0, [p, [EM H, S 
llll. [r3|. as well as in cold atomic gases [13], ID Kondo 
lattice 14 1 , and even nuclear and quark matter [3[ . The 
field has been energized by recent experimental evidence 
of the FFLO state in the heavy fermion superconduc- 
tor CeCoIn 5 [IH, [H, [ijj • Recent theoretical progress 
has indicated that the state may be achievable through 
applied current rather than an applied field 18]. The 
original prediction as well as the more recent ones rely 
on introducing a magnetic field (whether externally or 
as an effective internal mean field [l||), or using some 
other perturbation in order to create the unequal Fermi 
surfaces necessary to generate finite-momentum pairing. 

In this paper, a novel route to a 1D-FFLO state is 



proposed in the absence of a magnetic field or any other 
external perturbation with the AlAs system as an exam- 
ple. This system was recently fabricated and investigated 
by Moser et al. US HH ■ In such a system it should be 
possible to induce an FFLO state due to the interplay of 
Coulomb interaction effects with the specific bandstruc- 
ture of the AlAs quantum wire, which includes a unique 
inter-valley umklapp scattering process that is present at 
all densities. 
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FIG. 1: Schematic of the ID bandstructure in the AlAs quan- 
tum wire of Ref. I20I . The bandstructure has two degenerate 
subbands at the lowest densities, leading to four Fermi points. 
In addition, the band minima are separated by half an umk- 
lapp vector, ku/2, giving rise to umklapp interactions which 
are present at all fillings. The effective dispersion can be lin- 
earized so that the Fermi points, represented by black dots on 
the figure, are at k^ = -*£±k° F and k^ = ^±k° F where 
kjj is the umklapp vector, and k F is the magnitude of the 
Fermi wavevector measured from the bottom of each band. 
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A schematic of the AlAs quantum wire bandstructure 
is shown in Fig. [T] [2(1 SJ . While bulk AlAs has three de- 
generate bands, the growth direction and cleavage planes 
chosen in fabricating the wire leave two degenerate val- 
leys at low energy. In the effective one-dimensional band- 
structure, this leads to two degenerate nonoverlapping 
bands separated by half an umklapp vector, as shown in 
Fig. [TJ Because the momenta of the two band minima are 
related by half an umklapp vector, there is a class of umk- 
lapp excitations unique to this bandstructure which ex- 
ists at all densities. Because there are four Fermi points, 
the bandstructure supports several fermionic scattering 
processes. We retain only those scattering terms which 
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are relevant in a renormalization group sense. 




FIG. 2: Representative inter-valley umklapp scattering pro- 
cess which is responsible for encouraging the 1D-FFLO state. 
In the process which is shown, a pair of left and right movers 
from the A band scatters into a pair of left and right movers 
in the B band by absorbing an umklapp vector. This is an 
umklapp process which is present at all densities. 

In a quantum wire electrons are quantum mechanically 
confined to move along one direction with their motion 
in the remaining transverse directions confined via a po- 
tential, V r con /(rj_), where r± — (y,z) denotes the trans- 
verse coordinates of quantization. Electron-electron in- 
teractions within the wire are described by U (f) which is 
purely repulsive. The Hamiltonian is a sum of two inde- 
pendent terms in the transverse and longitudinal direc- 
tions with the result that the wavefunction (and there- 
fore the correlation functions) can be decomposed as a 
product of <j>(r±) and ip a (x) where (j>(f±) is the orthog- 
onal wavefunction of transverse quantization of the two 
degenerate bands (X and Y valleys) and ip s (x) the lon- 
gitudinal part. In order to describe the physics along the 
longitudinal direction we now promote the wavefunction, 
t/j s (x) , to the level of a field operator (for a field theoretic 
description) responsible for creating and annihilating the 
electrons taking part in the various scattering processes. 
With this in mind the second quantized Hamiltonian suit- 
able for the purposes of our study is 

H = J2f d 3 r^ s (r) ~ M + V con f(?x)) 

S ^ ' 

+ \YlJ d 3 rd 3 r'u(r~r')^l(r)^l,(r)^ sl (r )* s (r) 



(1) 



where ^ a (r) = ^(nJV'sfa) * s now ^ e field operator for 
an electron species of spin s — {t, 4}, and \i is the chem- 
ical potential in the leads. Because the low energy, long 
wavelength excitations occur in the vicinity of the Fermi 
points (see Fig. [2]) a further decomposition is possible 
with <J s (r) = <p(r ± )(ip As (x) + ip Bs (x)). The coordi- 
nate x is in the long direction of the wire. The longi- 
tudinal part of the field can be naturally expanded in 
terms of the right- and left- moving excitations, R ns (x) 
and L ns (x), respectively, residing around the Fermi 
points of the two bands (indicated by the black dots in 
Fig. [J) with ip As (x) = R A s{x)e lk r +x + L As (x)e lk r~ x and 
4>Bs{x) = R Bs (x)e ik P x + L Bs (x)e lk p x . The band in- 
dex is n — A, B and the Fermi momenta are defined by 



kp^ = —^f- ± kp and kp ± = Y ± k °F where k u is the 
umklapp vector, and k° F is the magnitude of the Fermi 
wavevector measured from the bottom of each band, as 
shown in Fig. [T] 

The low energy properties of the interacting ID elec- 
tron gas can now be conveniently described within the 
framework of the bosonization technique (22J. Within 
this approach, the fermionic field operators can be writ- 
ten in terms of bosonic fields <f) nv = (0 n f ± (j) n \)/y2 
and 9 n v = (#nt ± @nl ) / V2, where v = p,o (the charge 
and spin modes) correspond to the ± combination, s — 
{1,4} is the spin index and n — A, B is the band in- 
dex. The bosonic fields satisfy the commutation rela- 
tion [4>nv{x),d x , 6 n , v , {x )]=\n8{x-x )8 v y 8 n n , withfrset 
equal to one. We then have 

R ns ( x ) = VRns e yi[ e "p( g )-^( z )+ s ( fl ^( x )-'^(*))] 

v / 27ra 



(2) 



and, 



L n Jx) = VLns p % ^\- ^p( x )+^i'( x )+ s ( ^( x )+'t>^{x))] 
V27ra 

where a is the short distance cutoff, and ri Bns an d f]^ ns 
are the Klein factors for the right- and left- moving fields 
of band n with species of spin s. They are required to 
preserve the anti-commutation relations of the fermionic 
fields. The convenient field variables for the Hamiltonian 
in our problem will be a linear combination of the boson 
fields constructed out of the two bands. We define the 
transformation to a symmetric and an anti-symmetric 
basis as <f>± = -^{<f>Au ± §Bv) and 0± = -j%{0av ±0 Bl/ ). 

Upon bosonization and subsequent transformation the 
parts of the Hamiltonian corresponding to kinetic energy 
lead to harmonic terms in the symmetric and the anti- 
symmetric bosonic fields [<jr^ and 8^). The intraband in- 
teractions and the interband interactions, however, gen- 
erate both cosine interaction terms and harmonic terms 
[23l | . The Hamiltonian, H, can then be written in the 
following canonical form. 



H = 



IE 



U—p.rT 



dR 



v ±K±{d R e±Y+^{d R( j>t? 



Hi, 



(4) 

Here, R = (x + x )/2 is the center-of-mass coordinate of 
two electrons and a = x — x their relative coordinate in 
the long direction of the quantum wire. In the quadratic 
part of the Hamiltonian the bare symmetric and anti- 
symmetric Luttinger parameters, K„ can be expressed 
in terms of the original Luttinger parameters K u . The 
symmetric and anti-symmetric velocities v^j can also be 
expressed in terms of the original velocities v v . The terms 
responsible for interactions in the bosonic fields are in- 
cluded in Hi n t. There are eleven such interaction terms 
in the bosonized form. Weak coupling renormalization 
group (RG) (to second order in the interaction coupling 
constants) is then used to determine the instabilities of 
the AlAs quantum wire 0, [2f| H|| . The corresponding 
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coupled nonlinear differential equations are then solved, 
using appropriate physical inputs as initial conditions for 
the RG flow, as discussed below. 

The interaction terms which ultimately determine the 
final phases are 

t a 



100 



27r 2 a 2 

h 

27r 2 a 2 

to 

27r 2 a 2 
t d 



dR cos[v 47T0 ] cos[v 
dR cos[V47r(/)~] cos[V47r</)~] 
di?cos[\/47r6»p ] cos[\/47rV>+] 
dRcos[V^KB~] COs[a/47T^] 



27r 2 a 2 

The initial values of the coupling constants are set by 
t a = 2 J daM(a)cos[2k° F a] 
tb = 2 J daM(a) ( cos [2 k F a] — cos 
t c = 2 J daM(a) cos[2fc^a] cos 
t d = 2 M(a) cos 



(5) 



krra 



krra 



kjja 



(G) 



We use a screened Coulomb potential for the interac- 
tion kernel, M(a) = 4- — e , , , . . where e is the electronic 

' v I Aire y>w 2 +a 2 ' 

charge, e = 10.9e o is the dielectric constant of AlAs [27| . 
and e is the permittivity of free space. The width of the 
quantum wire w provides a short distance cutoff whereas 
the distance to the back-gate d provides a long distance 
cutoff. 

Using input parameters appropriate to the wire, and 
the expression for the coupling constants, Eq. [6j we can 
make an estimate for the initial conditions of the RG 
flow, i.e. the bare values of the cou pling co nstants U. 
The AlAs quantum wire of Moser et al. [20J, |2l| fabricated 
using standard design parameters has a transverse size of 
w w 15nm separated from the metallic gate by a distance 
d « 300nm. The wire length is L « 1/im, so that w/d « 
0.05 and L/d = 10/3. The parameters used to estimate 
k F for the AlAs bandstructure are as follows: the den- 
sity of the electrons in the quantum wire w 10 6 cm _1 and 
the effective mass of the electron m* — 0.3377i e along the 
long direction where m e is the bare mass of the electron 
[2l| . Experimental evidence [H, HI] suggests that spin 
rotational invariance is not broken in the AlAs quantum 
wells in the absence of a magnetic field (i.e. spin-orbit 
coupling is negligible.) Furthermore, since there is no ex- 
ternal magnetic field present spin rotational invariance is 
not broken and we set the initial spin Luttingcr parame- 
ters equal to one. 

Given these inputs for the bare parameters, we allow 
the running coupling constants to flow according to the 
RG equations with initial conditions for the bare Lut- 
tinger parameters K~ < 1 or K~ > 1 together with 
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FIG. 3: Variation of K~ with distance to the back gate (d) 
and width of the wire (w), using L — 1/ira, e = 10.9e o , and 
m* = 0.33m e . 



initial values for K~- set equal to one. When one or 
more coupling constants grows to be of 0(1), we stop 
the flow. At that point, the fields for the corresponding 
divergent couplings develop a gap to excitations, and ac- 
quire a definite expectation value. This expectation value 
is then inserted into the correlation functions in order to 
determine the type of phase this instability represents. 

Using this standard analysis, with initial conditions ap- 
propriate to the AlAs quantum wire as described above, 
we find that for the wire fabricated by Moser et al. the 
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coupling constants which diverge are t a and tb, leading to 
a state with divergent 2/c^-charge-density wave (CDW) 
correlations. This arises because the relative charge chan- 
nel flows to an effectively repulsive state, with effective 
Luttinger parameter K~* < 1. 

Although the bare interactions are (screened) repul- 
sive Coulomb interactions, it is possible for the interac- 
tions of the relative charge channel to flow to an effective 
attractive regime, i.e. K~* > 1. When this happens, 
the coupling constants which diverge are t c and t^, while 
t a and tb remain small. From the interaction terms we 
can then deduce that the dual antisymmetric charge field 
9~ gets gapped together with the spin fields <jr^. These 
gapped fields lead to a state with divergent intraband 
singlet superconductivity correlations with a finite pair- 
ing momentum - the 1D-FFLO state. The pairing 
mechanism in this case is shown in Fig. [2] The diverging 
interactions allow left and right movers from band "A" 
to scatter into a left and right mover in band "B" by 
absorbing an umklapp vector. That is, pairs form within 
each band in order to take advantage of the lowering of 
energy possible with umklapp scattering from one band 
to the other. This is analogous to the spin gap proximity 
effect mechanism described in Ref. |3Cj. Because each pair 
resides within a single valley, the pairs each have a net 
momentum, ±k v /2, leading to a ID FFLO state [Ull]. 

Phase diagrams showing this transition from 2k° F - 
CDW to the FFLO state are shown in Fig. [3 Our 



calculations indicate that the previously fabricated AlAs 
QWR of Refs. HJ would require extremely low den- 
sities (i.e. n as 10 4 cm _1 ) in order to achieve the FFLO 
state. Such low densities would not only make the behav- 
ior of the wire more susceptible to disorder, but it would 
also require a longer wire in order to achieve a reasonable 
total number of electrons in the wire. By using an alter- 
nate structure with a metallic side gate, it may be pos- 
sible to make the long distance screening length as small 
as d = 10 — 30 nm [33j. In this case, our calculations in- 
dicate that the FFLO state could be achievable for more 
reasonable densities, on the order of n « 10 5 cto _1 . 

In conclusion, we have proposed a novel route to an 
FFLO state in ID. The state should be achievable in 
AlAs QWR's for which the bandstructure leads to pair- 
ing generated by umklapp pair scattering which is present 
at all densities. This particular FFLO state is intrinsic 
to the wire, and does not require an external field or 
perturbation in order to induce the state. The theoreti- 
cal proposal is general. It is applicable to other systems 
with a similar bandstructure where analogous fermionic 
interaction terms are allowed. 

The author wishes to acknowledge Erica W. Carlson 
for suggesting the problem and additionally wishes to 
thank Matthew Grayson and Joel Moser for numerous 
insightful and useful discussions about the aluminum ar- 
senide quantum wire system. 
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